IFT-UAM/CSIC- 12-50; FTUAM-2012-92; HUPD-1201 



The string tension from smeared Wilson loops at large N 

Antonio Gonzalez-Arroyo a,b and Masanori Okawa' 
a Instituto de Fisica Teorica UAM/CSIC 
h Departamento de Fisica Teorica, C-15 
Universidad Autonoma de Madrid, E-28049-Madrid, Spain 
c Graduate School of Science, Hiroshima University, 
Higashi-Hiroshima, Hiroshima 739-8526, Japan £ 
(Dated: January 4, 2013) 

Abstract 

We present the results of a high statistics analysis of smeared Wilson loops in 4 dimensional SU(N) 
Yang-Mills theory for various values of N. The data is used to analyze the behaviour of smeared Creutz 
ratios, extracting from them the value of the string tension and other asymptotic parameters. A scaling 
analysis allows us to extrapolate to the continuum limit for N=3,5,6 and 8. The results are consistent with a 
1 IN 2 approach towards the large N limit. The same analysis is done for the TEK model (one-point lattice) 
for N=841 and a non-minimal symmetric twist with flux of k - 9. The results match perfectly with the 
extrapolated large N values, confirming the validity of the reduction idea for this range of parameters. 
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There is considerable interest in gauge theories at large N for their simplicity, proximity to phe- 
nomenologically interesting field theories and their presumed connection to string theory. Lattice 
gauge theory has proved to be a fundamental tool in deriving the non-perturbative properties of 
Yang-Mills theories at small N. In approaching large N, the standard pathway is to study the the- 
ory at increasing values of N and to extrapolate the results to infinite N. This is, no doubt, a costly 
procedure, with the additional risks involved in any extrapolation procedure. Nevertheless, results 
point towards a somewhat fast approach to the large N limit in many of its observables djLD — Q20 - An 
alternative would be to use the simplifications involved in the large N theory to find a way to simu- 
late it directly. An idea going in this direction is that of reduction or volume independenceiM-t4\- 



fl] 



[60. This allows the possibility of trading the space-time degrees of freedom with those of the 
group. The essential ingredient for the idea to work is invariance under Z 4 (N) symmetry, which is 
broken in the original proposal [3]-[4]. In the twisted Eguchi-Kawai model (TEK) [5], introduced 
by the present authors, an invariance subgroup is preserved at sufficiently weak coupling, enabling 
reduction to work. Recently, it was reported in Ref. [zl-[8]-[9] that symmetry-breaking takes 
place at intermediate couplings and N > 100. To circumvent this problem we proposed a slight 
variation of the model Q10Q. It exploits the freedom associated with an integer parameter entering 
the formulation, and representing the chromomagnetic flux through each two-dimensional plane. 
Traditionally this parameter was kept fixed when taking the large N limit, while we advocated the 
need to scale it with y/N in order to avoid symmetry-breaking phase transitions. In practice, the 
modification involves no additional technical or computational cost. Our initial tests [ | 1 0f | were free 
from the problems reported earlier. To further test the validity of this idea demanded performing 
state of the art computations of the large N observables and comparing them with those obtained 
for the TEK model. Furthermore, even if reduction operates at the level of the lattice model, our 
ultimate goal is the continuum theory, so a scaling analysis is necessary. These were our original 
motivations for embarking in the present work. 

Although other observables are possible, we have focused upon the string tension. This can be 
obtained as the slope of the linear quark-antiquark potential. Lately, the best determinations of the 
potential and of the string tension have been obtained by compactifying one dimension, and study- 
ing the connected correlation function of Polyakov lines D2D - In the large N limit this is subleading 
with respect to the disconnected term, and it is unclear how to make the connection. Thus, we stick 
to the traditional way in which the string tension is obtained from the expectation value of Wilson 
loops W(T,R). Here one meets a technical but severe difficulty, since large Wilson loops are very 



noisy quantities. Furthermore, the Wilson loops themselves are affected by ultraviolet divergences 
so that we will rather focus on the traditional Creutz ratios: 

,<rm i W(T + 0.5,R + 0.5)W(T-0.5,R-0.5) 

y(T,R) = -log (1) 

AK B W(T + 0.5,R-0.5)W(T -0.5, R + 0.5) 

which are defined for half-integer R and T. In the limit R « T these quantities are lattice approxi- 
mants to the force F(R) among quarks separated by a distance R. Although, Creutz ratios get rid of 
the constant and perimeter divergences in Wilson loops, they do so through a cancellation, which 
makes them even more numerically challenging. To reduce the errors we resort to the well-known 
Ape-smearing procedure[l 1] for the ordinary theory. The corresponding smearing for the TEK 
model is given by 

geared = Proj^ + c £( Z ^tf y t/„ t/J + l^U ,U v )\ (2) 

with z MV the twist tensor. Proj w stands for the operator that projects onto SU(N) matrices. This 
process can be iterated several times and produces a considerable noise reduction in the data. One 
could extract the string tension from the force F(R) obtained through Creutz ratios for R « T 
smeared in the three directions transverse to T. This is, however, very impractical in our case. 
It is much more effective to employ four-dimensional smearing and values R « T. The problem 
that arises in this approach is that not only the error, but also the value of the Wilson loops and 
Creutz ratios vary with the number of smearing steps. This could be an important source of 
systematic uncertainties, which might prevent a precision determination of the string tension from 
this source. To circumvent this problem, our strategy has been to use the smeared Creutz ratio 
values to extrapolate back and obtain un-smeared values. The extrapolated Creutz ratios do not 
depend on the number of smearing steps, and the errors are considerably smaller than the original 
un-smeared ratios. Having explained the main observables that we will be using, let us summarize 
in the next paragraph the goals and methodology used in this work. 

Our main goal is the determination of the string tension for large N Yang-Mills theory by 
means of the study of smeared Creutz ratios on the lattice. The large N value will be obtained 
by extrapolation of data taken at N=3,5,6,8 and by direct use of Twisted Eguchi-Kawai model at 
N = 841 and symmetric twist with flux k = 9. Indirectly, since the same procedures will be used 
to study the reduced and ordinary model, our results will serve to validate the reduction achieved 
by the TEK model in this physical range of parameters. Since the goal is the continuum result, we 
have simulated the model with Wilson action at several values of its coupling f3 = 2N 2 b = 2N 2 /A L , 
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FIG. 1. n s dependence of the Creutz ratio ^(4.5, 4.5) for N=841 and A L = 2.5974 

where A L is 't Hooft coupling on the lattice. The list of parameters and main lattice results are 
summarized in Table I. 

The analysis of data and presentation of results follows the following steps: 

1. Measurement of Wilson loops and Creutz ratios. 

For N=3, 5, 6, 8 lattice gauge theory, simulations are made on a 32 4 lattice with 260 configurations 
used for each A L . The number of configurations used in the TEK model for each A L is 5400, except 
at Ai=l.llll% where it is 2300. In both the ordinary theory and the TEK model, all configurations 
are separated by 100 sweeps, one sweep being defined by one-heat-bath update followed by five 
overrelaxation updates. 

We determined the Creutz ratios from Wilson loops smeared up to 20 times with c = 0. 1 in 
the range R,T e [3.5, 8.5]. Errors were estimated by jack-knife. Smaller values of R, T were also 
obtained, but dismissed for the analysis for being more sensitive to lattice artifacts. Larger loops 
can also be obtained but are too noisy and/or the number of smearing steps falls too short for them. 

2. Extrapolation to the un-smeared Creutz ratio with error. 

The extrapolation procedure depends on the values of R and T . For small values the smeared 
Creutz ratio are very well fitted to a dependence a(l - exp{-b/(n s + £)}) where n s is the number 
of smearing steps. This dependence is suggested by perturbation theory. In Fig. Q] we show, as an 
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0.63064 
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0.2453( 7) 
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0.63836 


0.01055( 12) 


0.2399( 14) 
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2.84625 


0.57441 


0.04028( 28) 


0.25 16( 40) 
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2.78676 


0.58892 


0.02668( 24) 


0.2501( 53) 
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2.76564 


0.59378 


0.02244( 20) 


0.2662( 15) 


5 


2.72242 


0.60338 


0.01654( 21) 


0.2700( 47) 
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2.65125 


0.61836 


0.01014( 10) 


0.2650( 15) 
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2.82408 


0.56997 


0.04126( 35) 


0.2580( 75) 
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0.58390 


0.027 11( 15) 
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0.2747( 29) 
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0.61363 


0.01090( 13) 


0.2657( 34) 
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2.80179 


0.56548 


0.04279( 26) 


0.2540( 19) 
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0.57930 


0.02840( 17) 


0.2672( 13) 
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2.72869 


0.58456 


0.02407( 14) 


0.2724( 13) 
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2.68902 


0.59405 


0.01811( 10) 


0.2728( 11) 
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2.62134 


0.60931 


0.01 112( 12) 


0.2726( 33) 


841 


2.77778 


0.55801 


0.04234(103) 


0.3019(170) 


841 


2.73973 


0.56902 


0.031 8 1( 60) 


0.2764(118) 


841 


2.70270 


0.57895 


0.02474( 56) 


0.2623(134) 


841 


2.66667 


0.58805 


0.01852( 45) 


0.2692( 94) 


841 


2.63158 


0.59651 


0.01418( 41) 


0.2722( 88) 


841 


2.59740 


0.60442 


0.01101( 24) 


0.2677( 49) 



TABLE I. We list the values of N and lattice couplings Al studied, together with the plaquette expectation 
value up and best fit parameters for Eq. |3] The N=841 case corresponds to the TEK model. 



example, the n s dependence of the R = T = 4.5 Creutz ratio for the TEK model and Al = 2.5914, 
together with the corresponding best fit curve. For larger values of R « T the first smearing 
steps represent a certain transient behaviour, which is then followed by a plateau, before decaying 
consistently with the previous formula. The extrapolated value is set to the plateau value. Details 
apart, it is important to emphasize that once the protocol to determine the un-smeared Creutz ratios 
was defined, it was applied by a program in exactly the same way for all values of N and Al, for 
both the reduced and ordinary model. Whenever bad fits or ambiguous behaviour was present, the 
errors were set to reflect the different options. 

3. Analysis of square R = T Creutz ratios 
The square Creutz ratios for large values of R = T are expected to behave as 

x (R,R) = K+ ^ + ... (3) 

A non-zero lattice string tension k is the consequence of Confinement. The linear term in 1 /R 2 is 
the predicted behaviour both from perturbation theory and from a string description of the quark- 
antiquark flux-tube. The dots contain corrections from different sources both of continuum and 
lattice origin. 

The data of both models and values of N show a very clear linear behaviour in 1 /R 2 , even at 
the smallest values of R. A good description of the data can be obtained with a three parameter fit 



based on Eq. [3j plus an additional term of the form tj/R 4 . The reduced chi square -\Jx 2 /ndf was 
typically of order 1 and never exceeded 2. The best fit parameters are listed in Table. I. 
4. Scaling analysis 

Since our goal is continuum physics we should extrapolate our results to the continuum limit. 
Scaling implies that, close enough to the continuum limit, results obtained at different values A L 
should coincide once the lattice spacing a(A L ) is chosen appropriately. In particular, the length of 
both sides of a rectangular Wilson loop are given in physical units by t = Ta(A L ) and r = Ra(AL). 
Using this fact and the definition of Creutz ratios one concludes: 

X (T,R) = a 2 (A L )F(t,r) + + + - ( 4 ) 

where the dots contain higher powers of a(A L ). The continuum function F(t, r) is given by 

. d 2 io g nyq,r) 

F{ur) = oVbT- (5) 

where 'Wit, r) is the value of the continuum t x r Wilson loop. Notice that, although the Wilson 
loop itself has perimeter and corner divergences, these disappear when taking the second derivative 



with respect to t and r. Thus, F(t, r) is a well-defined continuum function having the dimensions 
of energy square. 

In perturbation theory one gets 

*(f.r) = ypfe)(i + ±) (6) 

where jp is a given function of the aspect ratio z = r/t. For the full non-perturbative theory, one 
can study the behaviour of the function as t and r goes to infinity. One expects 

F(t,r) = cr + r(z)^ + ^j + ... (7) 

where cr is the string tension, and the dots represent subleading terms starting with l/(min(7, r)) 4 . 
The expansion is also exactly the same as predicted by an effective string theory description of the 
Wilson loop expectation value. 

This analysis justifies the parametrization used previously for square Creutz ratios with y = 
y(l) and k(A l ) = aa 2 {AL). In order to compute the continuum string tension we need to determine 
a(A L ). For very small values of At perturbation theory dictates its form: 

a(A L ) = -J- exp{-— 3— } (PoA^'V® = -^/W (8) 

However, it is well-known that scaling seems to work much beyond the region where Eq. |8] pro- 
vides a good approximation. There are several proposals in the literature, which have been dis- 
cussed and tested in many papers, which argue that Eq. [8] can be extended to the whole scaling 
region using improved couplings Ai(A£) in the previous formula, instead of A L itself. All these 
proposals can be considered perturbative renormalization prescriptions, and the ratio of the cor- 
responding scales is obtainable by a perturbative calculation, i.e. the ratio of lambda parameters. 

J~~L 

A particular proposal that has shown good results in previous studies was done by Parisi [|12] and 
used in the analysis of Ref. [13]. When expressed in A^ units it is given by a E = ^j-f(A E ), where 
Ae = (1 - up)SN 2 / (N 2 - 1). A somewhat different proposal resulted from the analysis of Allton 
et al. D2D. It is based on a different definition of the effective coupling Aa = Al/up (and a somewhat 
modified expression for f(A A )). 

Scaling then implies that the continuum string tension can be determined in A^ units as fol- 
lows: 




FIG. 2. N dependence of the continuum string tension 



Our data are consistent with the limit appearing in the right-hand side of the previous equation 
being approached linearly in a\. The extrapolated values for the ratio A^/ y/cr are displayed in 
Fig.[2]as a function of l/N 2 . Again a linear fit to the data with parameters 0.515(3) + 0.34(1)/W 2 
is quite satisfactory. The same procedure to obtain the string tension in the continuum limit was 
followed for the TEK model and N=841. The result, also displayed in the figure, is 0.513(6). The 
agreement with the large N extrapolated value of A^/ Vcr is very remarkable, and serves as a 
non-trivial test that reduction is operative for the TEK model in this range. 

Another remarkable feature of our result is that the N dependence matches perfectly with that 
obtained in Ref. (2D, which used different observables, techniques and range of 't Hooft couplings. 
The actual value of the large N ratio given in that reference was 0.503(2), which seems inconsistent 
with our result on statistical grounds. However, the estimated systematic errors quoted in Ref. fl 
are as large as 0.04. We should also mention that a recent analysis, largely complementary to ours, 
has obtained an estimate of the large N string tension which is consistent with our result [14]. 

In order to give a robust prediction for the large N string tension, we should also estimate 
our systematic errors. The most important source of these errors arises from an overall scale. If 
we repeat the procedure replacing the expression of by the formula given by Allton et al. B2D 
our estimate of the large N ratio A^/ ^/o- becomes 0.525(2). This is a 2 percent change in the 
predicted value, which is 5 times bigger than the statistical error. 

To give a more precise prediction we should use a non-perturbative renormalization prescription 
to fix a(A L ). It is possible to give a prescription based on Wilson loops and which follows the same 



philosophy as the one used to define the Sommer scale [15]. Let us consider the dimensionless 
function G(r) = r 2 F(r, r). A scale f can be defined as the one satisfying G(f) = G. If scaling holds, 
the choice of G is irrelevant (provided the equation has a solution), since it amounts to a change of 
units. For our analysis we took G = 1.65, by analogy with Sommer scale. However, we checked 
that taking other choices (G = 2. and G = 2.5) give consistent results up to a change of units. We 
recall that the idea of considering Creutz ratios with different aspect ratios z = R/T to define the 
scale appears in Ref. II 1611 . 

One possible way to determine the scale is by solving for R(A L ) in the equation 

G = R 2 (A L )x(R(A L ), R(A L ); A L ) (10) 

This gives us ci(Al) = f/R(A L ). Although, our data points are defined only for half integer R, it 
is easy to interpolate and obtain any real R. Interpolation is a much more robust procedure than 
extrapolation, and one can use different interpolating functions to estimate errors. 

The main problem of the previous procedure is that, as explained previously, the Creutz ratios 
have intrinsic scaling violations given by the second term in Eq. ©. Hence, a much better proce- 
dure is to make a simultaneous fit to all the square Creutz ratio data x(R, R) f° r a particular value 
of N (and all values of A£). Combining Eq. © and the expansion formula Eq. CO), one is led to the 
following functional form 

-2 -2 -4 2 

-MR,R\ A L ) = erf 2 + 2y r - + 4^(c + d%r) (11) 

Indeed, this formula describes remarkably well all our data, with chi squares per degree of free- 
dom of order 1 . Notice that the coefficient d accounts for the scaling violations arising from the 
definition of the Creutz ratios. The remaining terms parametrize the continuum function F(r, r) for 
large values of r. This depends on 3 parameters erf 2 , y and c. However, only two are independent 
since, by definition, f 2 F(f,f) = 1.65. This fixes Ac = 1.65 - erf 2 - 2y. If we substract the term 
proportional to d from the data, all data points should lie in a universal curve given by the function 
f 2 F(r, r) = ^G(r). In Fig. 3 we display the corresponding curve for SU(8) together with the best 
fit function extracted from Eq. (TTTT >. Notice how the values obtained from different couplings fall 
into a universal curve. Errors are displayed but hard to see at the scale of the graph. Similar curves 
are obtained for other values of Af and for the TEK model. 

The value of the parameters extracted from the fit are given in Table II. Notice that they are 
very similar for all theories. This makes the large N extrapolation very stable. A safe estimate of 
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TABLE II. The best fit parameters corresponding to Eq. QT] 



the large N value of erf 1 is 1.105(10), where the error now includes both statistical and systematic 
uncertainties. It is clear that a good part of the N dependence and systematic error found before 
resides in the ratio of scales rA^. 

In addition to the determination of the string tension, which sets the long-distance behaviour of 
Creutz ratios, there is considerable interest in the parameters that determine the approach to this 
long-distance limit. In particular, our results show that the large N slope parameter y has a value of 
0.272(5). The slope takes a non-zero value in perturbation theory equal to jp(\) = 2 (1 - 1 IN 1 ). 
Using this formula to define an effective coupling, our data implies ~ 8.4. At long distances, 
however, a new perspective arises which describes this term as arising from the fluctuation of the 
chromo-electric flux-tube stretching among the quark and the anti-quark. In the limit in which their 
separation is large compared to the thickness of this flux tube, an effective string theory description 
of the dynamics arises. The picture predicts (11711 that the coefficient of the 1 fr 1 contribution to the 
force F(r) is y(0) = f^. This prediction has been verified by lattice data. 

In our present case, it would be possible to study the function F(r, t) forri^t using information 
of non-square smeared Creutz ratios. In particular the function y(z) provides interesting infor- 
mation about the properties of the effective string theory. We can use our data to determine the 
function y(z) for the large N theory. The value for z = 1 coincides with the parameter y appearing 
in Table. II. Since by definition y{z) = y(l/z), we can parametrize this function in the vicinity of 
z = 1 as y(z) = y(l)(l + T^ip")- Our data for 7 > 0.5 allow a determination of r. For all values of 
N and A we get r = 0.31(6). 

As mentioned previously, the string picture predicts y(0) = n/12. However, the leading string 
fluctuation prediction for y(l) is ~ 0.16. Our numerical result for y(l) is far from this value and 
rather close to n/12. The same happens for the r coefficient, which is predicted to be close to 2. 
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FIG. 3. The continuum function r 2 F{r, r) = '—G{r) is plotted for SU(8) from our data. The solid lines 
correspond to the quadratic function Eq. QT] and to the linear part of this function. 



Remarkably, lowest order perturbation theory also has a prediction for r = 2/(n + 2) x 0.39, which 
is consistent with our data. The whole issue of string fluctuations for Wilson loops with different 
aspect ratios is being investigated at present [18]. 

In summary, we have presented a very precise measurement of the string tension for SU(N) 
Yang-Mills theory in the large N limit. It is remarkable that the N dependence is consistent with 
that obtained from correlation of Polyakov lines covering a different range of scales and distances 
r Vo" [2]. The large N result is also consistent with that obtained from the TEK single-site model, 
as predicted by the reduction idea. 
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